Chaotic dynamics of superconductor vortices in the plastic phase 



in 
o 
o 

(N 

o 

Q 



c ■ 
o ■ 

■ 

1 1 

5h ■ 
Or 

3 ; 



I 

C 

o 
o 



> 

o 
\o 
\o 
o 
in 
o 

-t— > 



i 

c 

o 
o 



dom pinning background in the (a;, y) plane. At T = 
the overdamped equation of motion of a vortex i in po- 
sition r, reads 



dt 



(^p) + F L (1) 
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We present numerical simulation results of driven vortex lattices in presence of random disorder 
at zero temperature. We show that the plastic dynamics is readily understood in the framework of 
chaos theory. Intermittency "routes to chaos" have been clearly identified, and positive Lyapunov 
exponents and broad-band noise, both characteristic of chaos, are found to coincide with the differ- 
ential resistance peak. Furthermore, the fractal dimension of the strange attractor reveals that the 
chaotic dynamics of vortices is low-dimensional. 

PACS numbers: 74.25.Qt, 74.40. +k, 05.45. Pq 

When flowing over a random medium, vortices in type 
II superconductors display a great variety of dynamical 
regimes, from the depinning threshold up to the high 
driving phase. Most of the V — I experiments 0, 0] and 
numerical simulations 0, 0, 0, IH reveal an intricate in- 
terplay between the 'peak effect' (PE), i.e. the increase of 
the depinning threshold current below the upper critical 
field H C 2, the peak of the differential resistance dV/dl, 
voltage noise and plastic flow of vortices. Below the PE 
an ordered phase is expected and the unusual excess noise 
measurements are understood within an edge contamina- 
tion process where a metastable disordered vortex phase 
generated at the edges is annealed into an ordered phase 
in the bulk Q . On the contrary, in the PE region a dis- 
ordered phase is expected and plasticity effects such as 
tearing are expected at the depinning threshold. These 
features have recently been studied in the mean field ap- 
proach 9] . However many open questions about the com- 
plex plastic flow remain, and in particular its dynamical 
properties. 

In this Letter, we propose to examine the plastic phase 
through the chaos theory of deterministic dissipative dy- 
namical systems. Charge density waves and Joseph- 
son junctions arrays have already been analysed through 
chaos theory 0], but such study is completely new for 
vortex lattices. We performed numerical simulations that 
clearly demonstrate the chaotic behavior of vortices in 
the plastic phase. While increasing the driving force, in- 
stabilities are developed by the non-linearities of the sys- 
tem and periodic regimes are destabilized giving rise to 
chaotic regimes. Such destabilizations have been clearly 
identified in our system to be the intermittency " route to 
chaos" . Furthermore, the broad-band voltage noise, the 
positive Lyapunov exponents and the fractal dimension of 
the strange attractor are used to characterize the chaotic 
phase which is shown to coincide with the peak of dV/ dl. 
A crucial result of our study shows that the chaotic dy- 
namics in the plastic phase is low-dimensional. Therefore 
within the framework of chaotic dynamical systems, our 
results open new perspectives in the understanding of 
vortex dynamics that are discussed in the conclusion. 

We consider N v Abrikosov vortices driven over a ran- 



where is the distance between vortices i and j, ri p 
is the distance between the vortex i and the pinning 
site located at r p , and V* is the 2D gradient opera- 
tor acting in the (x, y) plane. The viscosity coefficient 
is 77, F L = F L x is the Lorentz driving force due to 
an applied current. The vortex- vortex pairwise repul- 
sive interaction is given by a modified Bessel function 
U vv {rij) — 2eo A v Ko(rij / and the attractive pinning 
potential is given by U vp (ri P ) = —a p e~( Tip / Rp ) . In these 
expressions A v and a p are tunable parameters, Ax, is the 
magnetic penetration depth, and eo = (</>o/47rAx) 2 is an 
energy per unit length. We consider periodic boundary 
conditions of (L Xl L y ) sizes in the {x,y) plane. All de- 
tails about our method for computing the Bessel poten- 
tial with periodic conditions can be found in ^lj . Molec- 
ular Dynamics simulation is used for N v = 30 vortices 
in a rectangular basic cell (L x ,L y ) = (5, 6v3/2)Ax,- The 
number of pinning centers is set to N p = 30. We con- 
sider the London limit k = Ax/£ = 90, where £ is the 
superconducting coherence length |l2l | . The average vor- 
tex distance clq is set to ao = Ax, and R p = 0.22 Ax, 
T) = 1, A v = 2.83 x 10~ 3 Ax- We present results for two 
different pinning strengths corresponding to a maximum 
pinning force of F%> ax ~ 0.2F and F% ax ~ I.4F , where 
F = 2e A v /\L is a force defined by the Bessel interac- 
tion. In the weak pinning case, the driving force applied 
along a principal vortex lattice direction x is varied from 
up to F L - 3^0 ~ lOOi^ , where F C L is the critical 
Lorentz force along x. In the strong pinning case, the 
driving force is varied from up to F L ~ 3Fo ~ 20F C L . 
The choice of the double precision Runge-Kutta algo- 
rithm time iteration step St is dictated by the dominant 
force, and for example in the high driving phase we take 
St = 10-%, where t = r]\ L /F^ ax . 

The "experimental" procedure is the folllowing: we 
start by randomly throwing in the (x, y) plane N v vor- 
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Figure 1: Properties of the transition region from phase II 
to phase III in the weak pinning case. It clearly displays 
the type I intermittency route to chaos characteristics, a) 
part of the time evolution of the longitudinal velocity m (t) 
obtained for F L = 1.5221 x 1(T 4 . One sees laminar (i.e. 
periodic) phases interupted by chaotic bursts of large velocity 
fluctuations, b) distribution of the laminar phase durations of 
14 cm (t) measured for F L = 1.5221xl0" 4 . c) evolution of r max 
when varying the applied force F L from the intermittency 
threshold force F t . A clear power law is observed as shown 
with the line of slope —1/2. 

tices and N p gaussian pins, and relaxation with zero 
Lorentz force yields a vortex structure with dislocations. 
The Lorentz force is then slowly increased up to far in 
the high driving phase. The successive regimes we ob- 
serve in the weak pinning case F%f ax ~ 0.2F are the 
following. Phase I: pinned regime where all vortices have 
zero velocity. Phase II: plastic channels flowing through 
pinned regions and where the motion is either periodic or 
quasiperiodic as seen in Phase III: plastic flow with 
almost no stationnary vortices as seen in 0,01 El- In the 
following this motion shall be shown to be chaotic. Phase 
IV: fully elastic flow with no dislocations and where the 
motion occurs through rough static channels 
We shall now examine in details the plastic dynamics of 
vortices in the framework of classical chaos theory. The 
first central point of the letter shows that the transition 
from phase II to phase III is one of the three well-known 
"routes to chaos". In this very short applied force range, 
the typical longitudinal velocity of the vortex center of 
mass V x m that we measure in time is shown in Fig. la. 
It shows time intervals where the motion is periodic (the 
same as used to exist in phase II below the transition). 
The difference is now that such periodic regime becomes 
unstable and gives way to a chaotic burst displaying large 
velocity fluctuations. Then the system goes back to the 
periodic regime which is still unstable, giving way to an- 



other chaotic burst, and so on. The chaotic bursts corre- 
spond to apparently disordered trajectories of the moving 
vortices. However, from time to time, moving vortices are 
able to synchronize temporarily their motion into peri- 
odic motion (laminar phases). In the framework of the 
dissipative chaos theory, such intermittent regimes are 
known to be one possible way to drive the system from 
periodicity to chaos. The intermittency "route to chaos" 
has several characteristics and may be mainly classified 
in three types (I, II and III) depending on the unit cir- 
cle crossing value of the Floquet multipliers El 03. To 
determine the type of intermittency we observe in our 
system, we first measure for a given value of the applied 
force the distribution of the laminar (i.e. periodic) phase 
durations. Fig. lb shows such distribution obtained for 
V x m (t) displayed in Fig. la. This distribution of laminar 
phase durations shows a maximum at an upper bound 
i~max and a decrease for low durations which can be much 
smaller than T max . Furthermore, if we now increase very 
slowly the applied force in order to remain in an intermit- 
tent regime, the value of r max decreases as shown in Fig. 
lc. A very nice power law T max ~ (F L — F^)' 1 / 2 on al- 
most four decades is measured close to the intermittency 
threshold F t L , i.e. the force above which periodic regimes 
become unstable. The particular shape of the distribu- 
tion of the laminar phase durations and the exponent 
— 1/2 are characteristics of type I intermittency route to 
chaos related to a saddle-node bifurcation at F^ [TH , Hif . 
Note that the type of intermittency may change for dif- 
ferent pinning strengths (for stronger pinning parameters 
Fmax ~ 1-4-Fo, we observed for example a type II inter- 
mittency route to chaos characterized in particular by a 
different shape of the distribution of the laminar phase 
durations, and related to a subcritical Hopf bifurcation). 
Further increasing the applied force will give intermittent 
regimes with shorter laminar phase durations until they 
completly disappear therefore giving way only to large 
chaotic fluctuations. Then chaos expands in phase III. 

In the second central point of our letter, we examine 
in details the chaotic phase itself. Usual tools of chaos 
theory are successfully used to characterize the chaotic 
dynamics of vortices, and the link with the commonly 
observed differential resistance peak 0, 0> E| in vortex 
dynamics is established. This section characterizes the 
chaotic attractor of the vortices in the plastic phase III 
of the strong pinning case F^ ax ~ IAFq (where phases 
I, II and III are equivalent to those described above for 
the weak pinning case). We first compute the Lyapunov 
exponents. Positive Lyapunov exponents are a signa- 
ture of chaotic dynamics since they illustrate the "sen- 
sitive dependence on initial conditions" (SDIC) which 
is a property of chaotic attractors only. To compute 
the maximal Lyapunov exponent of our system we con- 
sider two very close initial conditions and observe how 
the distance d(t) in the phase space between the two 
corresponding trajectories evolve in time on the attrac- 
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Figure 2: Evolution of the maximal Lyapunov exponent A 
with the Lorentz force in phase III. Each point is the average 
of 20 couples of initial conditions and the error bars are the 
standard deviation. The inset displays the time evolution of 
the distance d(t) between two initial neighbouring trajectories 
in the phase space for F L = 0.0029 in the strong pinning case. 
The exponential divergence (positive slope A) characteristic of 
chaos is obvious for t < Tchaos- For time scales larger than 
Tchaos diffusive motion is observed as already shown in |g. 



tor. Since we integrate N v first order differential equa- 
tions of motion (Eq. QJ, the phase space is defined by 
the 2N V vortex coordinates and the distance d is defined 

'Xi(t) - Xdt)) 2 + (Yi(t) - ?i(t) 2 



by d 2 W = Eti 

where X t (t) = x l (t) - x cm {t), Y t (t) = y^t) - y cm {t), 
Xi(t) = Xi(t) - x cm (t), Yi(t) = yi(t) - y cm (t). In these 
expressions, (xi, y{) and (it, iji) are the vortex i coor- 
dinates, and (x cm , y cm ) and (x cm , y cm ) are the respec- 
tive coordinates of the center of mass. The tilde nota- 
tion (x,y) refers to the second trajectory generated by 
the neighbouring initial condition. The inset of Fig. [21 
displays an example of the time evolution of d we typ- 
ically find in phase III. It clearly shows an exponential 
divergence d ~ exp(Ai) of the two trajectories for time 
scales up to r c h a os ^ 1.6 x 10 4 ~ 50to- The slope there- 
fore defines a positive maximal Lyapunov exponent A 
characteristic of chaotic dynamics. Fig. [21 displays the 
evolution of A that we observe in phase III. The maxi- 
mal value of A therefore expresses that the fastest diver- 
gence of two chaotic trajectories occurs in the midrange 
of phase III. Finally, note that for time scales larger than 
Tchaos we find diffusive and superdiffusive motions (not 
shown) in the transverse and longitudinal directions as 
already reported in |(|. For a given applied force we 
now compute the power spectrum S a (f) of V£ m (t) , i.e. 
S a (f) = | III dt K m {t) exp(i2vft) | 2 where a = x 



Figure 3: Evolution of the strange attractor fractal dimension 
(computed with the correlation dimension v) with the Lorentz 
force in phase III. The inset displays the differential resistance 
curve (filled circles) and the low frequency longitudinal B x 
(open circles) and transverse B y (triangles) noises. Dotted 
lines separate phases II, III and IV. 



or y, and t% — t\ T c / laos . We define the low frequency 
noise B a by averaging S a {f) over the low frequency range 
. B x and B y are in some way a measure of the degree 
of chaos in the system since chaotic dynamics generates 
broad-band noise at low frequencies [16J . Concomitantly 
with the differential resistance dV£ m /dF L we plot in the 
inset of Fig. [3 the longitudinal B x and transverse B y low 
frequency noises. The rapid increasing of B x and B y con- 
firm the (rapid) setting of chaos in the vortex lattice and 
their maximal value coincide with the well-known peak 
of the differential resistance. As shown in Fig. [21 it also 
corresponds to the maximal value of A. Chaos is there- 
fore fully developed at the differential resistance peak. 
The positive Lyapunov exponents characrerize the ab- 
sence of temporal correlation in the chaotic regime due 
to SDIC. We shall now characterize the spatial corre- 
lations within the chaotic regime by computing the di- 
mension of the chaotic attractor. Such an attractor is 
known to be fractal with a non-integer dimension in the 
phase space. To characterize the fractal nature of this 
so-called strange attractor, we evaluate the correlation 
sum defined by C(p) = /im m ^ool/m 2 J]™ ;=1 i/(p - pti) 
which measures the number of couples of points [k, I) on 
the chaotic attractor whose distance pu is less than p. 
H(z) is the Heaviside function. For a limited range of 
p it is found that C(p) ~ p" where the exponent v is 
called the correlation dimension and is a simple measure 
of the fractal dimension of the attractor 0] . As shown in 
Fig- El the fractal dimension of the vortex strange attrac- 
tor has the same shape as the Lyapunov exponent curve 
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(Fig. [2J) and its maximum also coincides with the differ- 
ential resistance peak (inset of Fig.[3J). It therefore con- 
firms our previous finding that chaos is fully developed 
at the peak of the differential resistance curve. Above 
the peak, chaos still exists but its intensity decreases as 
shown by the deacreasing of v,\,B x and B y . It corre- 
sponds to the onset of transverse order of the smectic 
phase which we find to occur precisely at the peak. 
Furthermore, we find the very important result that the 
strange attractor of the driving vortices has a low (frac- 
tal) dimension, 1 < v < 2 (Fig.^J, which shows that the 
chaotic dynamics of a large number of vortices shrinks 
on a low-dimensional surface in the phase space. The 
crucial consequence of a strange attractor of dimension 
less than two is that the chaotic dynamics of vortices in 
the plastic phase may be described with only three dy- 
namical variables. This is an important result for further 
theoritical studies because simple analytical models with 
three dynamical variables should be sufficient to describe 
the complex plastic phase. Finally, we find that the end 
of chaos coincides with the dynamical freezing transition 
where the transverse velocity Vy" 1 drops to zero. We 
therefore clearly showed that the bottom of the differen- 
tial resistance peak marks the onset of chaos, while the 
end of the peak and the dynamical freezing transition 
appear as the end of chaos. 

In conclusion, we obtained conclusive results about the 
chaotic dynamics of vortices in the plastic phase. The 
route to chaos has been identified in details: type I (II) 
intermittency in the weak (strong) pinning case is found. 
Chaos characterized by positive Lyapunov exponents and 
broad band noise is found to coincide with the differen- 
tial resistance peak. Furthermore, the fractal dimension 
of the strange attractor shows that the chaotic dynamics 
of vortices is low-dimensional. Therefore our results open 
new perspectives in the theoretical understanding of the 
plastic flow phase which is much less developed than the 
fast moving vortex phases 0, 0] and than the plastic 
dcpinning transition |9( . In particular we show the impor- 
tant result that the chaotic dynamics in the plastic phase 
may be understood with only three dynamical variables. 
Hence, our results combined with the usual tools of the 
chaos theory (e.g. bifurcation theory, embedding dimen- 
sions, Poincare sections, strange attractors, time series 
analysis) should help for further theoretical, numerical 
and experimental investigations of the open issues re- 
lated to the formation of fractal objects in the complex 
space phase of driven plastic vortices. Finally, our results 
let foresee new possibilities of controlling vortex motion 
for device applications using the concept of controlling 
chaos developed these last years (see e.g. 0). The goal 
of such control procedure is to lock the chaotic system 
into a stable periodic orbit which either used to be un- 
stable in the uncontrolled system (feedback scheme by 
weakly changing parameters |2l|). or which is newly cre- 
ated (non-feedback scheme by weakly forcing the system, 



see e.g. |22j). We therefore suggest that the concept of 
controlling chaos might be used to design device applica- 
tions to rectify plastic (chaotic) vortex motion. 
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